Geometry of the energy landscape of the self-gravitating ring 
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We study the global geometry of the energy landscape of a simple model of a self-gravitating 
system, the self-gravitating ring (SGR). This is done by endowing the configuration space with a 
metric such that the dynamical trajectories are identified with geodesies. The average curvature 
and curvature fluctuations of the energy landscape are computed by means of Monte Carlo simu- 
^ 1 lations and, when possible, of a mean-field method, showing that these global geometric quantities 

X^/y . provide a clear geometric characterization of the collapse phase transition occurring in the SGR 

as the transition from a flat landscape at high energies to a landscape with mainly positive but 
Ov fluctuating curvature in the collapsed phase. Moreover, curvature fluctuations show a maximum in 

correspondence with the energy of a possible further transition, occurring at lower energies than the 
' '. collapse one, whose existence had been previously conjectured on the basis of a local analysis of the 

energy landscape and whose effect on the usual thermodynamic quantities, if any, is extremely weak. 
We also estimate the largest Lyapunov exponent A of the SGR using the geometric observables. The 
geometric estimate always gives the correct order of magnitude of A and is also quantitatively correct 
at small energy densities and, in the limit N — >• oo, in the whole homogeneous phase. 



PACS numbers: 05.20.-y; 05.70.Fh; 02.40.-k; 05.45.-a; 95.10.Fh 



I. INTRODUCTION 

Given a system of classical interacting particles described by the Hamiltonian 
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all the information on the dynamics as well as on the equilibrium collective properties of the system is encoded in 
the potential energy V. Indeed, the forces entering the equations of motion are given by the gradient of V, and the 
equilibrium statistical measures in both the canonical ad the microcanonical ensembles can be defined in terms of V, 
since the kinetic energy contribution can be integrated out explicitly. 

To extract information from the function V one may use a class of methods commonly referred to as "energy 
landscape" methods [l[. The value of the function V(q) at a given configuration q — (gi, . . . ,qjv) € M, where M 
is the iV-dimensional configuration space, is the height of the landscape, so that configurations with low energy 
are valleys (and their bottoms are minima of V), maxima of V correspond to peaks, an saddle points of V are 
mountain passes[5l|. Hence the topography of the energy landscape is determined by the stationary points of V, 
i.e. the configurations q s such that W(q s ) — 0. Starting from the pioneering work by Stillinger and Weber 0, the 
energy landscape approach has been fruitfully applied to study glassy systems (see e.g. Q and references therein) and 
biologically motivated problems such as protein folding (see e.g. Q and references therein). The basic idea behind the 
energy landscape 52] methods is very simple, yet powerful: if a system has a rugged, complex energy landscape, with 
many minima and valleys separated by barriers of different height, its dynamics will experience a variety of time scales, 
with oscillations in the valleys and jumps from one valley to another. One of the drawbacks of this approach is that 
one should in principle find all the stationary points of V to obtain a complete characterization of the landscape, or 
at least all the minima and all the saddles connecting them: for many-particle systems and generic potential energies 
this is an impossible task. Also assuming that a partial sampling of the landscape is sufficient, a huge computational 
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effort is required in order to find a non-negiigible fraction of minima and saddles. It is then natural to ask whether 
a characterization of the energy landscape can be obtained in terms of global quantities, averaged on the landscape 
itself (see 0-13 for a related discussion in the context of toy models of proteinlike polymers). 

The energy landscape approach, with its emphasis on the topography of the basins around minima, has a topological 
flavour [53j. However, the potential energy V can also induce a geometric structure on the configuration space: a 
metric function g on M can be defined in terms of the potential energy V such that the configuration space becomes a 
Riemannian, or pseudo-Riemannian, manifold. We can thus speak of the geometry of the energy landscape, and as we 
shall see in the following the averages of geometric quantities can provide the global characterization of the landscape 
we are looking for. The intuitive reason why geometric information on the landscape, and especially curvature, could 
be a relevant one is that the dynamics on a landscape would be heavily affected by the local curvature: minima of the 
energy landscape are associated to positive curvatures and stable dynamics, while saddles involve negative curvatures, 
at least along some direction, thus implying some instability. One can reasonably expect that the arrangement and 
detailed properties of minima and saddles might reflect in some global feature of the distribution of curvatures of the 
landscape, when averaged along a typical trajectory. 

There are many possible ways to define a metric structure on the energy landscape. The most immediate choice 
would probably be that of considering as our manifold M not the configuration space but the TV-dimensional surface 
z = V(qi, . . . , <7iv) itself, i.e., the graph of the potential energy V as a function of the N coordinates qi,...,qn of 
the configuration space, and to define g as the metric induced on that surface by its immersion in R JV + 1 , Although 
perfectly reasonable, this choice has two major drawbacks. First, the explicit expressions for the geometric quantities 
in terms of V and its derivatives are rather complicated. Second, and most relevant to what we are going to discuss 
in the following, with that choice of the metric the link between the properties of the dynamics and the geometry 
is not very precise, i.e., one cannot prove that the geometry completely determines the dynamics and its stability. 
We shall thus choose a metric g, referred to as the Eisenhart metric such that the dynamical trajectories are the 
geodesies of the configuration space endowed with such metric and that the stability of the dynamics is completely 
determined by the curvature of such metric. In addition, the explicit expressions for the curvature of the Eisenhart 
metric in terms of derivatives of V are particularly simple (see Sec. |TT]). As reviewed in [9j-fllj|. with such metric one 
can define a suitable curvature observable whose average and fluctuation allow to characterize the global properties 
of the energy landscape as well as to estimate the degree of instability of the dynamics, quantified by the largest 
Lyapunov exponent. 

The aim of this paper is to apply these geometric tools to characterize the energy landscape and the dynamics of 
a simple model of a self-gravitating system, the self-gravitating ring (SGR), introduced in [12|]. The microcanonical 
statistical mechanics of this model has been studied by a mean- field technique in [ljj], showing that the system 
undergoes a phase transition between a collapsed and a uniform phase. The latter phase transition has been related 
to some particular stationary points of the energy landscape in [14| , where also some indication of the possibility of 
another phase transition, occurring at lower energies, was presented. In the following, we shall compute the curvature 
of the energy landscape defined via the Eisenhart metric and show that its average and fluctuations clearly mark the 
collapse phase transitions; moreover, the fluctuations of the curvature also display a peak at the energy where the 
second phase transition conjectured in [14| should occur. We shall also estimate the largest Lyapunov exponent in 
terms of the average and fluctuations of the curvature and compare such estimates with direct numerical calculations. 
The geometric quantities provide a reasonable estimate of the Lyapunov exponents for any energy, and an accurate 
one for large and small energies; at intermediate energies the geometric estimate, though still giving the correct order 
of magnitude, is worse and we shall discuss the origin of this result. 

The paper is organized as follows: in Sec. [TT] we shall briefly recall the relation between geometry and dynamics, 
in order to define the geometric quantities that will be studied to characterize the global properties of the landscape 
and that will be used to estimate the largest Lyapunov exponent. In Sec. IHII we shall describe the model studied, i.e., 
the self-gravitating ring (SGR), and its equilibrium statistical properties; Sec. II VI will be devoted to presenting and 
discussing our results on the geometric quantities as well as on the Lyapunov exponents, and we shall end the paper 
with some concluding remarks in Sec. [V] 

II. GEOMETRY AND DYNAMICS 

That the dynamical trajectories of a classical dynamical system whose Hamiltonian is given by Eq. (JTJ) can be viewed 
as geodesies of the configuration space endowed with a suitable Riemannian metric is a classic result in analytical 
mechanics (see e.g. [15l|). The first to conjecture that this correspondence may be relevant to the physics of many-body 
systems — and especially to the foundations of statistical mechanics — was probably N. S. Krylov, who in his doctoral 
thesis of 1942 (reprinted in [lj|) suggested that the dynamical instability of the trajectories (what we nowadays call 
Hamiltonian chaos) may be due to the negative curvatures of the configuration manifold. Krylov's suggestion was 
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based on earlier mathematical results 1 1 7L 1 1 811 ; since then, geodesies flows on negatively curved manifolds have become 
an important chapter of ergodic theory [19[ . However, in many dynamical systems of physical interest curvatures do 
not have a definite sign or even are only positive, yet the dynamics appears to be chaotic. In these cases a subtler 
mechanism may be at work: the oscillations of the curvature alonga dynamical trajectories may induce chaos in close 
analogy to parametric instability, as first suggested by Pettini [20| . Starting from this observation the existence of 
different chaotic regimes in many Hamiltonian systems was related to the curvature of the configuration space 0, HH 
and an analytical estimate of the largest Lyap unov exponent, that is the commonly used measure of the strength of 
chaos in a dynamical system, was derived [22l |23| . 

In the following we shall very briefly recall the main definitions and results necessary to understand the application 
of the geometric approach to the SGR model: in Sec. Ill Al we shall describe the use of the Eisenhart metric and the 
derivation of an effective stability equation allowing to estimate the Lyapunov exponent, while in Sec. Ill Bl we shall 
discuss how such an equation describes an effective geometry of the energy landscape. Full details and references to 
related work on different models can be found in the reviews 13-Hll. 



A. Eisenhart metric and curvature 



As already argued above and as it will become clear in the following, an important geometric quantity to characterize 
the global properties of the energy landscape is curvature. The definition of the curvature of a manifold M depends 
on the choice of a metric g [3, Hfl : once the couple (M, g) is given, a covariant derivative V and a curvature tensor 



R(e 



can be defined; the latter measures the noncommutativity of the covariant derivatives in the coordinate 



directions e$ and ej . A scalar measure of the curvature at any given point q £ M is the the sectional curvature 

K(ei,ej) = (i?(e i ,e i )e i ,e i ) , 



(2) 



where (•, •) stands for the scalar product. At any point of an iV-dimensional manifold there are N(N — 1) sectional 
curvatures, whose knowledge determines the full curvature tensor at that point. One can however define some simpler 
curvatures (paying the price of losing some information): the Ricci curvature Kn(ei) is the sum of the if's over the 
N — 1 directions orthogonal to e^, 



N 

K R {ei) =22K(ei,ej) , 
and summing also on the N directions e, one gets the scalar curvature 

N N 
i=l i,j=l 

then, -j^j and N ^_^ can be considered as average curvatures at a given point 



(3) 
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As argued in previous works (see e.g. [|| and references therein), the Eisenhart metric @ turns out to be a 
particularly good choice. Given a system with Hamiltonian (JT|) and setting m = 1 without loss of generality, this 
(pseudo-Riemannian) metric is defined on a configuration space with two extra dimensions, M x R 2 , with local 
coordinates ((fa, Qi, <ZiV) <?jv+i)) an d its arc-length is 



ds 2 = S i , j dq i dq j - 2V(q){dq ) 2 + 2dq°dq N+1 
The metric tensor will be referred to as gE and its components are 
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as can be derived from Eq. (|5|). 

The geodesies, i.e. the "straight lines" on the manifold, are defined by the equation 



V„u = , 



(7) 
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where v is the velocity vector along the geodesies itself. This means that geodesies are the curves whose (covariant) 
acceleration vanishes. Writing Eq. (J7J in local coordinates with the metric gE one finds that the geodesies coincide 
with the dynamical trajectories of Hamiltonian (JlJ, provided one restricts to the subset of geodesies such that the 
parametrization of their arc-length is affine, 

ds 2 = C 2 dt 2 , (8) 
where we can set C 2 = 1 without loss of generality. The nonvanishing components of the curvature tensor are 

d 2 V 

where i, j — 1, . . . , N. It can then be shown that the Ricci curvature ([3]) in the direction of motion, i.e., in the direction 
of the velocity vector v of the geodesic, is given by 

k r {v) = Ay , (io) 

where AV is the Laplacian of the potential V, and that the scalar curvature 1Z identically vanishes. We note 
that Kr(v) is nothing but a scalar measure of the average curvature "felt" by the system during its evolution; we 
will refer to it simply as Kn dropping the dependence on the direction. Another feature of Kr is its very simple 
analytical expression which simplifies both analytical calculation and numerical estimates. It is also worth noticing 
that expression (fTQ| is a very natural and intuitive measure of the curvature of the energy landscape, as it can be 
seen as a naive generalization of the curvature f"(x) of the graph of a one- variable function to the graph of the 
TV-dimensional function V(qi, ?jv)' the Laplacian of the function. However, the previous discussion shows that it 
is much more than a naive measure of curvature and that it contains information on the local neighborhood of the 
dynamical trajectories, i.e., on their stability. As we shall see in the following, this information is sufficient to extract 
an estimate of the largest Lyapunov exponent of the system. 



1. Geometric estimate of the Lyapunov exponent 

The stability of the geodesies ([7]) is completely determined by the curvature of the manifold: given a geodesic whose 
velocity is v, a vector field J, referred to as the Jacobi field, can be defined such that (J, v) = and it measures the 
distance between nearby geodesies. The Jacobi field obeys the Jacobi equation [25( 

VlJ = R(v,J)v. (11) 

Written in local coordinates, Eq. (| 1 1 [) becomes (summation over repeated indices is understood) 

™+#*J>% (12) 
ds 2 JM ds ds y J 

and in the particular case of the Eisenhart metric we get @ 

d 2 J* d 2 V Jk 

dP + dq*dq k ' [ ' 

which is commonly referred to as the tangent dynamics equation: the largest Lyapunov exponent is the mean expo- 
nential growth rate of the solutions of the latter equation. Hence using the Eisenhart metric the growth rate of the 
Jacobi field is directly related to the Lyapunov exponent of the system. 

In general Eq. (|lip is very complicated. There are anyhow two special cases where it becomes simple. The first 
case is when the sectional curvature ([2]) of the manifold is a constant; such manifolds arc referred to as isotropic 
manifolds [54,]. Choosing for simplicity a geodesic reference frame, that is a frame parallelly transported along the 
geodesic, the Jacobi equation can be written 



d 2 J l 
ds 2 



KJ l = 0, (14) 



where J 1 is any of the components of J in the chosen frame, s is the arclength and K is the constant sectional 
curvature. The character of the solutions depends on the sign of K . If K > 0, J remains bounded and oscillates 
with frequency \f~K. If K = 0, J grows linearly. Hence on a positive (or zero) constant curvature manifold all the 
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geodesies are stable (or marginally stable). On the contrary, on a manifold with constant negative curvature all the 
geodesies are unstable, since the solutions of Eq. ([T^| with K < grow exponentially, J(s) oc exp (y/—Ks) . All the 
eodesics are unstable also on a manifold with non-constant sectional curvature, provided the latter is strictly negative 
25} , and this explains why geodesic hows on negatively curved manifolds play such a relevant role in ergodic theory 
19]; unfortunately, no similarly general result can be proved for manifolds of non constant positive curvature (or for 
manifolds whose curvature does not have a definite sign). The other case where Eq. (|11|) becomes simple enough is 
that of two-dimensional manifolds, without any constraints on their curvature. In this case the Jacobi equation in a 
geodesic frame reads as 

^L + K(s)J = 0, (15) 

where J is the single component of the Jacobi field (remember that J is orthogonal to v so that the geodesic frame 
can be chosen such as J has components only along one direction) and K(s) = K(v, J) is the sectional curvature in 
the directions of v and J. We see that Eq. (|15[) is remarkably similar to the Jacobi equation for constant-curvature 
manifolds, but for that K(s) can now vary along the geodesic. The solutions of Eq. (TT5t can grow exponentially, 
implying instability of the geodesies, for two reasons: either K < always, or the oscillations of K induce, via a 
sort of parametric instability, an exponential growth of the envelope of J(s). As first suggested in (2pj . where also 
some numerical examples were worked out, and further elaborated in [2ll423| . this mechanism can be at work also in 
high-dimensional manifolds and can explain the origin of chaos in those systems whose curvatures are mainly positive. 

Starting from this idea, an effective scalar stability equation can be derived from the Jacobi equation (ITT1) under 
the assumption of quasi-isotropy of the manifold @ . Loosely speaking, quasi- isotropy means that the oscillations of 
the curvature along a geodesic are small as compared to the average curvature along the geodesic itself. Without 
entering technical details, when N —> oo this assumption and the use of the Eisenhart metric allow to approximate 
Eq. pT|) with the following scalar equation 



k{t)ip = 0, (16) 



where t is time[55], ip denotes any of the components of an effective Jacobi field, and k(t) is an effective curvature 
(defined below). Equation (|16p is formally very close to Eq. (|15p. with a remarkable difference: it is a stochastic 
equation. The effective curvature k(t) is a Gaussian stationary stochastic process whose probability distribution 
is defined by the equilibrium distribution of the Ricci curvature Kr per degree of freedom, as calculated in the 
microcanonical ensemble. More precisely, the average of k(t) is given by 

*o = (k) = ^(K R )^ = i(AV)p (17) 

and its fluctuation is 

at = (k 2 ) (fc) 2 = 1 ({AV 2 )» - (AV)l) . (18) 

In Eqs. (|17p and (ITS)) , stands for the microcanonical average. Moreover, kit) is a <5-correlated process: 

(fc(ti)fc(ta)) - (k{h)){k{t 2 )) = Ta 2 k S(t 2 - tO , (19) 

where the correlation time scale r can be self-consistently estimated as follows @. First, two independent time scales 
7i and r% are identified as, respectively, 



and 



n = 9/ _- (20) 

2yfco + ok 



r 2 = ^. (21) 
07c 



The scale T\ is the dominant time scale as long as <C k ; for — > it is related to the period of the oscillations 
of the Jacobi field on a manifold with constant (positive) curvature fco- The other scale T2 dominates when the 
fluctuations grow, eventually until they become of the same order of magnitude of the average curvature. Then r is 
obtained from t\ and Ti such that the shortest time scale dominates: 



T- l =T^+T^. (22) 
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Having completely defined the process k(t), the average exponential growth rate of the solutions of Eq. (fTB]) can 
be calculated [1^] analytically [H, [23|. Since, as recalled above, the Jacobi equation written for the Eisenhart metric 
coincides with the tangent dynamics equation (fl3|) . this directly yields a geometric estimate of the largest Lyapunov 
exponent A as 



A(fc ,CTfc,T) 



A 



2 V 3A 




1/3 



(23a) 
(23b) 



Some remarks are in order as to Eq. (|23p. First, it yields an estimate of the largest Lyapunov exponent in terms of the 
average and fluctuations of the effective curvature, once r is defined as a function of fco and Ofc via Eq. (|22[) . without 
any free parameter. However, the above definition of r is by no means a direct consequence of any theoretical result, 
but only a rough, physically based estimate. We have thus explicitly indicated the dependence of A on r in Eq. (|2"3")l 
because such an estimate might well be improved independently of the general geometric framework. In [26| it was 
also considered as a free parameter to be fitted after comparison with numerical data. In this work we shall keep the 
estimate (|22[) for t in order to have a parameter-free prediction for the Lyapunov exponent A. Second, the ingredients 
fco and entering Eq. (|23p can be calculated analytically whenever the microcanonical average and fluctuations of 
the Laplacian of the potential energy of the system under investigation can be calculated exactly; otherwise, they 
can be obtained via relatively short numerical simulations, in any case much shorter than the long integrations of 
the tangent dynamics equation (fT5]> needed to directly calculate Lyapunov exponents, especially when the latter are 
small. Finally, when 07. <C ko, that is when the quasi-isotropy assumption is more likely to be correct, Eq. (|2"3"|) yields 



A 



' k 5 



(24) 



thus elucidating the deep relation between curvature fluctuations and dynamical instability. 

In many cases Eq. (|2"3"|) yields good, if not very good, estimates of the largest Lyapunov exponent [9j; in the case of 
the Fermi-Pasta-Ulam j3- model (FPU-/?) , that is a chain of oscillators with quartic nonlinearity, it yields a prediction 
whose agreement with numerical data is excellent over eight orders of magnitude in energy and six orders of magnitude 
in A [l,|23j]. Being the theory based on the quasi-isotropy assumption, one expects it may fail when Ok > fco- Indeed, 
the quantitative agreement between the estimate (|23[) and numerically calculated Lyapunov exponents is typically 
worse when the fluctuations are large as compared to the average, although the order of magnitude of the estimate is 
often correct as well as the estimated dependence of A on the energy density [H, ii. 



B. Curvature and collective behaviour: the effective geometry of the landscape 

Apart from giving the basis of the Lyapunov exponent estimate ([2"3"| . the effective stability equation ([TST lsuggests 
that the dynamics of the system can be seen as a geodesic flow on an effective two-dimensional manifold [57] whose 
curvature fluctuates around an average value fco , the amplitude of the curvature fluctuations being . This provides 
an intuitive geometric picture of the energy landscape: when the curvature fluctuations are small with respect to the 
average curvature, the landscape is very smooth and one expects the dynamics to be only weakly chaotic, while if 
the fluctuations of the curvature are of he same order of magnitude as the average the landscape is rough and chaos 
should be strong. Indeed, this is what happens in quite a few Hamiltonian systems with many degrees of freedom, the 
FPU-/? model being the paradigmatic example: a crossover between a weakly and a strongly chaotic regime, which 
has been referred to as the strong stochasticity threshold (SST) [27, 28], occurs when the energy density e = E/N is 
increased and it does corres pon d to a crossover between a smooth effective manifold, where Uk <C fco, and a rougher 
one, with cr fe /fco ~ 0(1) [Hip. 

Even more dramatic changes in the effective geometry of the landscape have been detected in systems undergoing 
thermodynamic phase transitions: for instance, in models with short-range interactions where the phase transition 
breaks a 0(n) symmetry, curvature fluctuations <7fc exhibit a cusp- like behaviour at the transition @, [l(J |2fj [29l43Tj| . 
This suggests that a geometric counterpart of the phase transition may be found in the dependence of the effective 
geometry of the energy landscape on the energy density. A case of particular interest is that of the ferromagnetic mean- 
field XY model, often referred to as the hamiltonian mean- field (HMF) model [Hj], that is a model of fully coupled 
planar classical spins with ferromagnetic interactions, where a mean-field ferromagnetic-paramagnetic transition occurs 
while increasing the energy density above a critical value e c . In this model the quantities fco and Cfe describing the 
effective geometry of the landscape can be calculated analytically, as first done by Firpo in [33|] to estimate the 
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Lyapunov exponent using Eq. (|23|) (see also Ref. Q for a discussion of that result from the point of view of the 
present discussion). The calculation shows that the phase transition corresponds to a geometric transition from 
a manifold with positive curvature and nonvanishing curvature fluctuations to a flat manifold with zero curvature 
fluctuations. The effective average curvature fco decreases with energy and vanishes continuously at the transition, 
while the curvature fluctuations increase with e as long as e < e c and discontinuously jump to zero at e c . Both 
ko and au vanish for e > e c . As we shall see in the following, this model coincides with a particular limit of the SGR 
model we are going to study, so that it is of particular interest to the present work. 

Before describing the SGR model and studying the effective geometry of its energy landscape, it is worth briefly 
mentioning some further developments which have been triggered by the geometric approach. First, constructing 
some abstract geometric models able to reproduce the cusp-like behaviour of the curvature fluctuations at a phase 
transition [29j and noticing that such behaviours occur also with different metric functions [30j . it was conjectured 
that these geometric transitions were a consequence of deeper topology changes: this was explicity shown in the 
particular case of the mean-field XY model [HJ and quite a few research activity followed (see [ll|, |35[ and references 
therein). Second, the study of the effective energy landscape geometry of a toy model of a polymer has shown that 
the behaviour of curvature fluctuations allows to discriminate between polymers that have a proteinlike behaviour, 
collapsing to a well-defined native state, and polymers that collapse to a random state 0-0|- Hence this tool seems 
to be useful also to characterize transitional phenomena in intrinsically finite systems like proteinlike polymers. 



III. THE MODEL 



In this Section we shall describe the self-gravitating ring (SGR) model and its equilibrium statistical properties, in 
particular in the microcanonical ensemble. In Sec. IIII A"l we shall discuss the model; in Sec. IIII 51 we shall discuss its 
equilibrium statistical behavior. 



A. The self-gravitating ring 

The SGR model describes N identical classical particles of mass m constrained to move on a ring of radius R and 
mutually interacting via a regularized gravitational potential. It was introduced in [l2j and its Hamiltonian is 

i N i N ^2 

1 v 2 1 v Gm 

= h 1 " 2 U ec* (*-*,)+ a ' (25) 

where i?j £ (— 7r,7r] is the angular coordinate of the i-th particle, 

P l= mR 2 -^ 26 
at 

is its angular momentum, R is the radius of the circle and G is Newton's gravitational constant. The dimensionless 
constant a > provides the short-distance regularization: when i?j -1? 3 <1 the interaction between the i-th and the 
j-th particle is effectively harmonic on a length scale of the order of 

d a = RV2a . (27) 

The SGR model can thus describe a self- gravitating system when d a <C R, i.e., a -C 1. On the other hand, in the 
opposite limit a — > oo the SGR model becomes equivalent [5 8] to the ferromagnetic mean- field XY model, or HMF 
model. 

Following [TH [3(| we can define a characteristic time scale r (the typical time a particle needs to go around the 
circle when all the others are uniformly distributed) as 



R 3 
GNm 

and dimensionless momenta as 



(28) 



Vi = ~^, i = l,...,N. (29) 
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One can thus write a dimensionless Hamiltonian % as 

r 2 1 N 1 N 1 

H = ^ SGR -2^~^-r 2 £ ^-cos^-^ + a " (30) 

The dimensionless Hamiltonian Q30p is extensive, due to the jj rescaling of the coupling between the particles. The 
Kac prescription for making extensive a long-range interaction [37| is here obtained via a suitable choice of the 
adimensionalization procedure. From now on we will consider only the dimensionless Hamiltonian (|30[) . 



B. Microcanonical thermodynamics 

The equilibrium thermodynamics of the SGR has been studied in 13] and later in [38[ by a mean-field method, 
expected to be exact in the thermodynamic limit N — > oo. The main features of the thermodynamics depend on 
the choice of the softening parameter a and can be summarized as follows. There is a phase transition between a 
homogeneous high-energy phase and an inhomogeneous low-energy phase, occuring at a critical energy density e c 
weakly dependent on a (for instance, e c ~ —0.227 for a = 1 and e c ~ —0.32 for a — 10~ 2 ). In the microcanonical 
ensemble the phase transition is discontinuous if a < 10 -4 and continuous otherwise, i.e., there is a microcanonical 
tricritical point at a ~ 10~ 4 . There is inequivalence between the canonical and the microcanonical ensemble (and 
the system exhibits a region of microcanonical negative specific heat) for any a < 10 _1 ; for larger values of a, the 
transition is continuous in both ensembles, the region of microcanonical negative specific heat disappears, there is 
ensemble equivalence and the thermodynamics of the SGR is qualitatively similar to that of the mean-field XY model. 

In the rest of the paper we shall consider three values of a: a = 1, where the thermodynamics is very close to that 
of the mean-field XY model; a = 10 -2 , where there is a region of negative specific heat but the microcanonical phase 
transition is still continuous; and finally a — 3 x where the microcanonical behaviour of the SGR is closer to that 
of a real self-gravitating system, with a large region of negative specific heat and a discontinuous phase transition. 
This choice of values of a allows to explore all the qualitatively different regimes exhibited by the svstem; moreover, 
for these values results on the thermodynamicSjUseful for validating our results, were obtained in [38[. 

Since the mean-field technique developed in [13| does not allow to calculate fluctuations of thermodynamic observ- 
ables, to c omp ute Cfc we had to resort to Monte Carlo (MC) simulations, performed using the Ray microcanonical 
algorithm [391 ] . In Figs. 11121 and 131 we plot the caloric curves of the SGR as obtained with the mean- field method in 
[38j for the three chosen values of the softening parameter a and we compare them with the results obtained using 
the MC algorithm with N ranging from 50 to 400, showing a very good agreement and the absence of appreciable N- 
dependence of the MC results in the inhomogeneous phase e < e c . We note that MC results extend to smaller energies 
than those explored with the mean-field technique: the latter indeed does not easily converge for small energies. 



IV. RESULTS AND DISCUSSION 



Let us now present and discuss the results obtained for the average effective curvature ko and curvature fluctuations 
<7fc of the energy landscape of the SGR. We shall first use these quantities to give a geometric characterization of the 
phase transition and then to estimate the largest Lyapunov exponent. Both quantities have been obtained by means 
of the MC method already mentioned in Sec. IIII Bl the average curvature ko has been calculated also using the 
mean-field technique of Ref. [l3j |. yielding results valid in the N — > oo limit. The number of particles considered in 
the MC simulations, as in the case of the computation of the caloric curves, ranged from N = 50 to N = 400. The 
SGR Hamiltonian (|30|) is long-ranged and we did not use any cutoff or approximation on the forces as those typically 
used to simulate self-gravitating systems in astrophysical contexts (see e.g. [4(| SJ), so that the computational cost 
grows as iV 2 and this sets a limit to the maximum number of particles that can be simulated in reasonable time. The 
computation of the caloric curves has shown that at least in the inhomogeneous phase the iV-dependence of the MC 
results is practically negligible, so that we may expect that these numbers of particles are sufficient to yield a good 
estimate of the geometric properties of the landscape at large AT. In all simulation runs we used O(10 6 ) MC steps per 
particle. 

In Sec. IIV Al we present the results for ko and <Jk and discuss how they provide a geometric interpretation of 
the homogeneous-inhomogeneous phase transition. In Sec. IIV A II we argue that these results may also support the 
existence of another transition in the SGR, at lower energies than s c , as suggested in [TJ]. Finally, in Sec. IIV Bl we 
use ko and er^ to estimate the largest Lyapunov exponent of the SGR and we compare the theoretical estimates with 
direct numerical calculations. 



9 

0.25 

0.2 
0.15 

0.1 
0.05 


-0.35 -0.3 -0.25 e c -0.2 -0.15 

e 

FIG. 1: (Color online) Microcanonical caloric curve T(e) of the SGR model for a = 1. The blue solid line has been obtained 
in [3^] using the mean-field technique developed in [l3|; the symbols are MC results for N = 50 (green circles), N = 100 (cyan 
rhombs), N = 200 (red squares), and N = 400 (blue triangles). Statistical errors are smaller than the symbols' size. The 
dotted vertical line marks the critical energy density. 
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FIG. 2: (Color online) As in Fig.[T]for a = 10" 2 . 

A. Geometry and phase transitions 

The average effective curvature fco and fluctuations ak for the SGR are given by inserting the potential energy of 
the Hamiltonian into the definitions (|17l) and (TT51) ; as already mentioned, the microcanonical averages have been 
obtained by MC simulations (and also using the mean- field method for fco, at sufficiently high energies). The results 
for these quantities are reported in Figs. 21 [5j and[U for a = 1, a = 10~ 2 , and a = 3x 10 -5 , respectively. 

Starting from the case a = 1 (Fig. we observe that the average curvature fco is positive, decreases with the 
energy density e and vanishes continuously at e c ; it remains equal to zero in the whole homogeneous phase. More 
precisely, the mean- field method, valid as N — > 00, yields fco = Ve > e c ; the MC data for e > e c are always positive 
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FIG. 3: (Color online) As in Fig. U for a = 3 x 10 
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but very small and exhibit a clear tendence to vanish as N grows. For e < e c the MC data agree with the mean-field 
solution as long as the latter one can be reliably computed, and do not display relevant finite-size effects (in any 
case, these effects are smaller than in the homogeneous phase) as in the case of the caloric curves. The behaviour of 
the curvature fluctuations <ik is more complicated: grows with e until a maximum is reached at e max below e c , 
then it starts decreasing smoothly until, close to e c , a sharp decrease sets in and for s > e c the data show a strong 
iV-dependence, compatible with a% = for N — > oo; indeed, in the homogeneous phase we have at(N) cx N~ s with 
5 ~ 0.5. The data thus suggest that, as N — ¥ oo, <Jk jumps discontinuously from a finite value to zero at e c . The 
qualitative behaviour of kg and ak is thus very close to that observed in the mean-field XY model (see Sec. IIV Al) . 
the only relevant difference being that in the latter ak reaches its maximum right at e c where it also jumps to zero, 
while in the SGR the maximum occurs at e m ax < Ec We have anyway considered also larger values of a (data not 
shown) and observed that e max moves towards e c as a grows, thus suggesting that at large values of a, and ideally 
when a — > oo, the behaviour of ak is exactly that of the mean-field XY model, as expected. 

Considering now the two smaller values of a (Figs. [5] and [5]), as far as fco is considered the qualitative behaviour 
is essentially the same as before, apart from a slight change in the shape of the curve k (e): k > for e < e c and 
k$ — for £ > e c . However, while the transition between the two regimes seems continuous when a = 10 -2 (Fig. 
[5]), our data are consistent with a discontinuous jump of ko at s c for the smallest value of a we considered (Fig. [6]). 
The behaviour of the the curvature fluctuations 07. is instead very different from the a = 1 case. Here, after reaching 
a maximum at an energy density £ max well below e c (e m ax of the same order of magnitude of, but not equal to, the 
energy density where the specific heat changes sign from positive to negative) <7fc(e) smoothly decreases to almost 
zero, and the residual jump at s c is either absent, or, if present, is negligible as compared to the maximum value 
cr(e max ). In the homogeneous phase, also for these values of a we have ak{N) oc N~ s with S ~ 0.5. Apart from the 
already mentioned possible jump of ko at the transition for a = 3 x 10 -5 , no other appreciable differences between 
a = 10~ 2 and a = 3 x 10~ 5 show up in the graphs of fco(e) and afe(e). There is a big change in the numerical values, 
that is however to be ascribed to the changes of the energy and length scales induced by the variation of a. Hence 
a clear difference in the effective curvature of the energy landscape emerges between those values of a where the 
specific heat is always positive and there is equivalence between the canonical and microcanonical ensemble and those 
values where a negative specific heat region exists and the canonical and microcanonical descriptions are no longer 
equivalent. Moreover, the fact that the microcanonical phase transition is discontinuous seems to be reflected in a 
discontinuity of kg at e c when a = 3x 10~ 5 . 

The global geometric picture of the energy landscape that emerges from these results is very similar to that of 
the mean-field XY model: the phase transition from an inhomogeneous to a homogeneous state corresponds to a 
transition from a positive-curvature manifold, with energy-dependent curvature fluctuations, to a flat manifold with 
no curvature fluctuations. The manifold corresponding to the homogeneous phase is easily identified with the whole 
configuration space, which is an A^-torus T^, endowed with an effectively flat metric. This is true regardless of the 
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FIG. 4: (Color online) Average effective curvature ko (top panel) and effective curvature fluctuation <jfc (bottom panel) of the 
energy landscape of the SGR as a function of the energy density e for a — 1. Symbols refer to MC results for TV — 50 (green 
circles), TV = 100 (cyan rhombs), TV = 200 (red squares), and TV = 400 (blue triangles). Errorbars indicate statistical errors. 
The solid blue line in the upper panel is the mean-field result for TV — > oo. The dotted vertical line marks the critical energy 
density. 



value of a. The differences emerge in the collapsed phase, and the data we have obtained suggest that these differences 
are related to how this torus is "constructed" starting from an energy landscape whose geometry is that of a sphere 
(positive curvature with vanishing curvature fluctuations) at small energy densities. Starting from the minimum 
energy density, when the latter increases the average curvature k decreases while the curvature fluctuations increase 
until they reach a maximum at £ max . The fact that fco decreases while Ofc increases may be related to that larger 
and larger regions with smaller — or even zero or negative — local curvature are visited by the system and the effective 
manifold becomes somewhat intermediate between a sphere and a torus. This happens for any value of a. Then, 
after the maximum the fluctuations starts to decrease; if a is sufficiently large, they suddenly jump to zero, as if the 
rest of the torus were constructed all of a sudden, by pinching all the remaining holes all together. Indeed, in Ref. 



34] a topological analysis of the configuration space of the mean-field XY model showed that in the inhomogeneous 



phase a "half-torus" is progressively built, which becomes a full torus all of a sudden at the transition. For smaller 
values of a this sharp geometric transition is no longer present (or at least makes the curvature fluctuations change 
only by a small fraction of their maximum value): also the rest of the torus is progressively built, as its first half. The 
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FIG. 5: (Color online) As in Fig. H for a = 1CT 2 . 



presence of a negative specific heat in thermodynamics appears thus to be related to a region where the "second half" 
of the torus, with more and more regions of zero or even negative curvatures, is progressively explored: this region 
is "compressed" in a sharper geometric change when the specific heat is always positive. It must be noted that if 
ascribing the decrease of the average curvature to the visiting of larger and larger regions of zero or negative curvature 
is correct, one would naturally expect that in the homogeneous phase the systems visits both positive and negative 
curvatures: if that occurs with equal probability, fco vanishes but cr^ reaches a constant value, instead of vanishing. 
Indeed, this is what happens in XY models with short range- interactions [HI, |3l| , whose configuration space is also 
T N . In our case we have fco = but also = 0, so that it seems that in the homogeneous phase the system visits 
only regions with zero curvatures, or that the metric becomes fiat on the whole manifold. The fact that the effective 
energy landscape becomes flat in the homogeneous phase is thus an effect of the long-range nature of the interactions. 

As a final remark, it is important to note that for all the values of a there exists an energy region in the inho- 
mogeneous phase where au > fco- The latter region extends from an energy density slightly smaller than e max to e c , 
as shown in Fig. where the extrapolation of the results on fco and <7fc to N — > oo is reported. The presence of 
such a region has two main consequences. First, it means that the effective geometry of the landscape is particularly 
complicated in that region, since fluctuations of the curvature larger than the average curvature itself show up on the 
macroscopic scale. Second, in this region the hypothesis of quasi-isotropy of the manifold at the basis of the geometric 
estimate (|23|) of the Lyapunov exponent A no longer holds, and we may expect this reflects in a disagreement between 
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FIG. 6: (Color online) As in Fig. [4] for a — 3 x 10 5 . The inset in the top panel shows a magnified view close to the transition 
to appreciate the discontinuity in ko at e c suggested by the mean- field results (blue solid line). 



estimated and directly calculated values of A. As we shall see in Sec. lIVBl this is precisely what happens. But before 
considering Lyapunov exponents, let us concentrate on another aspect of the effective geometry of the landscape. 



1. Another transition? 



In Ref. (14| the stationary points of the energy landscape of the SGR have been studied, showing that the 
homogeneous-inhomogeneous phase transition is related to the fact that a class of stationary points becomes asymp- 
totically flat, that is the determinant of the Hessian of the potential energy vanishes for those stationary points when 
N — > oo. The asymptotic flatness had been proposed as a necessity criterion for selecting stationary points able to 
induce a phase transition in [42|, [43[ and has been later on referred to as the KSS criterion. Being based on stationary 
points of the potential energy landscape, the KSS criterion singles out a value of the potential energy density v — V/N 
where a phase transition can occur. In Ref. flij it was shown that not only v — 0, that is the maximum of the potential 
energy density, corresponding to the homogeneous phase, satisfies the KSS criterion, but also another value, which 
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FIG. 7: (Color online) Estimate of the effective geometry of the energy landscape of the SGR extrapolated to N — » oo as a 
function of the energy density e for a = 1 (top panel), a = 10 -2 (middle panel), and a = 3 x 10~ 5 (bottom panel). The blue 
solid line is the average effective curvature ko and the red dashed line is the effective curvature fluctuation at- The curves have 
been obtained by mean-field results where available (that is, for ko for all the values of e but the smaller ones), by interpolating 
the data with N — 400 in the inhomogeneous phase and by setting at to zero in the homogeneous phase. The dotted vertical 
line marks the critical energy density. 
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FIG. 8: (Color online) Comparison between the potential energy density at the maximum of the curvature fluctuations u ma x 
defined in Eq. (|32[) (red symbols) and the potential energy density v c predicted by the KSS criterion as given by Eq. (|31[) (blue 
solid line) as a function of the softening parameter a. Inset: the same in a log-log scale, to show a wider range of values of a. 



depends on the softening parameter a: 

4 + a[6 + a(5 + 2a)] 



2a [(2 + a) 3 / 2 + a 3 / 2 ] 



(31) 



Hence another phase transition might be present in the SGR, at a critical energy such that the average value of 
the potential energy density equals v c . It was also suggested that this transition, if present, could separate a low- 
energy region where the probability that a particle goes around the circle drops to zero (so that the system has a 
core-only density profile) from a higher-energy region where the probability of exploring the whole circle becomes 
nonzero, although the phase is still clustered (so that the system has a core-halo density profile). Such a change in the 
distribution may have very small effects on the thermodynamics. In a recent work [44j no signs of such a transition 
have been detected in the thermodynamic observables of the SGR, but indeed a change in the density profile of the 
same kind as anticipated above has been observed. A similar behaviour has been observed in a toy model that well 
approximates the SGR when a is small (see Ref. [36| , where it was also shown that the same happens in the Thirring 
model g3). 

We are now going to show that at v c something happens also in the effective geometry of the energy landscape, and 
precisely a maximum in the curvature fluctuations. We have already noted that the curvature fluctuations show a 
maximum at an energy density £ max which is always smaller than e c but tends to the latter value when a — > oo. On 
the other hand, previous studies have shown that a maximum in the curvature fluctuations appears to be often (if 
not always) associated to big changes in the collective behaviour of the system: typically to a thermodynamic phase 
transition [9j (in this case the curvature fluctuations appear to be singular at the maximum, e.g. showing a a cusp) 
or to a big configurational change which is not a thermodynamic transition due to the finiteness of the system, as in 
the toy model of a proteinlike polymer studied in 0, 0] . In the latter case the maximum of the curvature fluctuations 
is smooth. The maximum of the curvature fluctuations of the SGR does not show a tendency to become singular as 
TV grows, so that in a sense it seems more similar to that of the polymer. It remains to show that the maximum of 
the curvature fluctuations occurs at v c \ in Fig. [5] we plot the microcanonical average of the potential energy density 
at the energy density £ ma x of the maximum of the curvature fluctuations, 

"max = Jf( V )» (£max) , (32) 

as a function of the softening parameter a, together with the predicted value v c given by Eq. (|3"Tj) . The agreement 
is quite good. This result suggests that a deep change in the properties of the system, which as observed in (44| 
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FIG. 9: (Color online) Geometric estimate of the largest Lyapunov exponent A of the SGR as a function of the energy density 
e for a = 1, obtained inserting MC results for fco and into Eq. (|23p with N = 50 (green circles), N = 100 (cyan rhombs), 
TV = 200 (red squares), and N = 400 (blue triangles). Errorbars indicate statistical errors. The dotted vertical line marks the 
critical energy density. 



is probably not a usual thermodynamic transition but a change from a core-only to a core-halo density profile, 
actually occurs at v c . We also observe that this confirms once again the sensitivity of the effective geometry of the 
energy landscape to changes in the collective properties which may not affect usual thermodynamic observables like 
temperature, internal energy or specific heat, as it happens for proteinlike polymers [1, H|. 



B. Geometry and chaos 



We can now use the results on the average effective curvature fco and curvature fluctuations presented in Sec. 
IIV Al to estimate the largest Lyapunov exponent of the SGR using Eq. (|23]l. 

First of all, we note that our results, obtained with N < 400, are consistent with the expectation that the dynamics 
is regular in the thermodynamic limit in the homogeneous phase, i.e., A — > when N — s- oo and e > e c , since both 
fco and A vanish in the thermodynamic limit. More precisely, using Eq. (j23p. when e > e c we obtained X(N) cx TV -7 , 
with 7 ~ 0.28 for a = 1, 7 ~ 0.31 for a = 10 -2 , and 7 ~ 0.38 for a = 3 x 10 -5 . These iV-dependencies compare 
reasonably well with theoretical estimate made in the case of the mean- field XY model (33|, that is 7 = 1/3. At 
variance with the homogeneous phase, the iV-dependence of the geometric estimate of A in the inhomogeneous phase 
e < e c is weak, as a consequence of the weak iV-dependence of fc and o>. We can thus assume that the estimate of A 
obtained inserting in Eq. ([2"B")) our MC results for fco and Ofc with N = 400 is a reasonable estimate also for N —> 00, 
in the collapsed phase. In Fig.|H]we show the geometric estimate of A obtained from Eq. (|2"3"|) for a = 1 and N ranging 
from 50 to 400. 

In Figs. fT0l[TT1 and[T2lwe compare the geometric estimates of A for a = 1, a = 10~ 2 and a = 3x 10~ 5 , respectively, 
with direct measurements of A obtained by means of molecular dynamics (MD) simulations. For the sake of clarity, 
only the results for the largest number of particles considered are reported for each value of a, since the geometric 
theory is expected to perform better for larger systems [Hj]. The MD simulations were performed by numerically 
integrating both the equations of motion derived from the SGR Hamiltonian (|30[) and the tangent dynamics equation 
(|13p . i.e., the equation of motion in the tangent space obtained by linearizing the equations of motion al ong the 
dynamical trajectory. Numerical integration was performed using the McLachlan-Atela symplectic integrator [46j and 
relative fluctuations of the energy remained smaller than 10 -5 in the worst case. The computational cost of the MD 
simulations of the SGR grows as iV 2 as the MC one, but is even higher because to compute Lyapunov exponents one 
has to simulate both the dynamics in the configuration space and that in the tangent space. Moreover to extract 
reliable values of A with the standard method [47| very long simulations are needed. We were thus able to obtain 
reliable direct numerical estimates of A for systems up to N = 400 with a = 1 and a = 10~ 2 , but only up to N = 100 
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FIG. 10: (Color online) Comparison between the geometric estimate (|23[) of the largest Lyapunov exponent A (blue squares) 
and the direct measurement of A obtained by molecular dynamics simulations (red circles), as a function of the energy density 
e for a — 1 and N = 400. Errorbars indicate statistical errors. The dotted vertical line marks the critical energy density. 



with a = 3 x 1CP 5 . Although only the results for the largest iV's are shown, also the direct MD estimates of A are 
consistent with A — > when N —> oo and s > e c , at least for the two larger values of a considered; for a = 1 we found 
X(N) oc TV -0 - 2 and X(N) oc iV~ - 25 for a = 10~ 2 . For a = 3x 10 -5 we were unable to extract the TV-dependence of 
A. We also note that a residual TV-dependence of the MD values of A shows up also for e < e c , at variance with the 
geometric estimates. The latter TV-dependence is larger at smaller a's. For small e, no appreciable TV-dependence was 
observed. 

Several observations are in order as Figs. [TOl [TT] and [12] are concerned. First, the order of magnitude of the 
geometric estimate (|23|) always compares very well with the MD data, for any a and any e, the worst difference 
between geometric estimates and MD data being a factor of 1.5 for a = 1 and a factor of 2 for a = 1CP 2 and 
a = 3 x 10~ 5 . This is far from being trivial, since the absolute values of A vary by more than a factor of 100 
passing from a = 1 to a = 3 x 10~ 5 . The same can be said of the overall behaviour of the curves A(e), which is 
reasonably reproduced by the geometric estimates, although the agreement is worse as a gets smaller. Second, the 
agreement between geometric estimates and MD data becomes quantitatively very good at small e for all the a's 
considered, and also in the whole homogeneous phase e > s c if we assume that the extrapolation A — > for large 
systems is correct. These two situations are precisely those where the theory leading to Eq. (|23l) is expected to be 
correct. In the homogeneous phase the correctness of the geometric theory as TV — > oo is somewhat trivial: geodesies 
of a flat manifold are not chaotic. However, also at finite TV the geometric estimate yields reasonable results in the 
homogeneous phase. At small e one has <Jk <C ko so that the quasi-isotropy assumption is consistent and the theory 
performs very well as expected. Third, in an intermediate energy region there is a quantitative disagreement between 
theory and direct numerical measurements. As already pointed out in Sec. IIV Al in this region > ko so that the 
quasi-isotropy assumption is questionable and a disagreement is not surprising at all, and has been observed in all 
the cases where curvature fluctuations grow too much with respect to the average curvature (see e.g. Q). The more 
so, because <Jk > ko means that negative curvatures appear with non-negligible probability, and the combined effect 
of the parametric instability due to fluctuations and the defocusing effect due to negative curvature is often very 
difficult to understand, resulting sometimes in a reinforcement of the overall chaoticity and sometimes in an inhibition 
of chaos (see e.g. [Hj]). In some cases as that of the 1-d XY model [23[ the geometric estimate of A can be corrected 
to account for the negative curvatures' contribution, yielding a very good agreement between theory and numerics 
also when ak > fco- We did not try to implement a correction as such in the present case, nor did we try to use the 
correlation timescale t, estimated in Eq. (|2"2")h as a fitting parameter as done in [26| . 

Finally, we note that the quantitative disagreement between theoretical estimates and numerical measurements of 
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FIG. 11: (Color online) As in Fig. \M with a = 10" 




FIG. 12: (Color online) As in Fig. [TO] with a = 3 x 10" 5 and N = 100. 



A when a > kp may also have a different origin. As correctly pointed out by Vallejos and Anteneodo in [48| and 
also by Politi |49|, Eq. (|23|) does not directly estimate the Lyapunov exponent A but rather a generalized Lyapunov 
exponent commonly referred to as A2. This is due to that, for ease of calculation, in deriving Eq. (|23|) one first 
averages the trajectories of the stochastic oscillator over the realizations of the stochastic process and then computes 
the exponential growth rate of the averaged solution. To properly compute A one should instead first compute the 
exponential growth rate of the solutions and then take the average over the realizations (this is formally very similar 
to the difference between quenched and annealed averages in the statistical physics of disordered systems). When the 
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amplitude of the fluctuations is small A2 is typically a very good estimate of A, but this is not always the case when 
the amplitude of the fluctuations get larger. However, the pursuit of this idea is left for future work. 



V. SUMMARY AND CONCLUDING REMARKS 



We have studied the global geometry of the energy landscape of the SGR by endowing the configuration space with 
the Eisenhart metric, such that the dynamical trajectories are identified with geodesies. The average curvature and 
curvature fluctuations of the energy landscape have been computed by means of MC simulations and, when possible, 
of a mean-field method, showing that these global geometric quantities provide a geometric characterization of the 
collapse phase transition. In particular, the average curvature ko as a function of the energy density e behaves as 
an order parameter, vanishing for e > e c . Curvature fluctuations au provide further insight into the nature of the 
collapsed phase and on its dependence on the softening parameter a. The behaviour of Ufc(e) also shows a maximum in 
correspondence with the energy of a possible further transition, occurring at lower energies, whose existence had been 
previously conjectured on the basis of a local analysis of the energy landscape [l4|. Such a transition may correspond to 
a change from a a core-only to a core-halo density profile, although it does not appear to have appreciable consequences 
on the usual thermodynamic quantities This confirms once more the sensitivity of the geometric observables to 
changes in the collective properties which may not affect usual thermodynamic observables like temperature, internal 
energy or specific heat, as it happens for instance in the case of the proteinlike polymers considered in [HQ. 

The effective curvature and its fluctuations allowed us to estimate the largest Lyapunov exponent A of the SGR 
via Eq. ([23]) . The geometric estimate always gives the correct order of magnitude of A and is also quantitatively 
correct at small e and, in the limit N — > 00, in the whole homogeneous phase e > e c . For intermediate energies 
the accuracy of the estimate is worse. Though not surprising at all, since for these values of the energy density the 
assumptions leading to Eq. (|23p may fail, this result nonetheless calls for a deeper understanding which may also 
lead to an improvement of the accuracy of the geometric approach to chaos, as discussed in Sec. IIVB1 Even if, in our 
opinion, the main value of the geometric approach is in its explanation of the deep origin of chaos in Hamiltonian 
many-particle systems and not in its effectiveness as a computational tool for Lyapunov exponents, it is also true that 
any theory is valid inasmuch as it can predict the outcomes of (numerical, in the present case) experiments. Hence 
the possibility of modifying Eq. (l23l) to incorporate a better estimate of A is a very interesting suggestion for future 
work. 

To conclude, all the results on the effective geometry of the energy landscape of the SGR show how this global 
approach to the properties of the energy landscape may be useful and complementary to the local analysis of stationary 
points, as that carried out in [TJ] for the same model. 
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